


Institutional Archive of the Naval Postgraduate School 


Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1971 


Adaptive angle tracking and correlation for 
airborne direction-finding. 


Pfendtner, Felix. 


Monterey, California ; Naval Postgraduate School 


http://ndl.handle.net/10945/15555 


Downloaded from NPS Archive: Calhoun 


| Calhoun is the Naval Postgraduate School's public access digital repository for 
D U DLEY research materials and institutional publications created by the NPS community. 
sa Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS'‘s first 
KNOX appointed — and published — scholarly author. 





LIBRARY Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 


http://www.nps.edu/library Monterey, California USA 93943 





— ————— = a 


ADAPTIVE ANGLE TRACKING AND CORRELATION 
FOR AIRBORNE DIRECTION-FINDING 


by 


Felix Pfendtner 

















United States | 
Naval Postgraduate Schoc! 


by 


FOR AIRBORNE DIRECTION-FINDING 
Felix Pfendtner 





Thesis Advisor Poe hd Eas 
March 1971 
eH 


Approved for pupolic release; distrubution unlimcted. 





Adaptive Angle Trackine and Correlation 


for Airborne Direction-Finding 


by 


Felix Pfendtner 
Lieutenant Commander, Federal German Navy 


Submitted in partial fulfillment of the. 
requirements for the degree of 


ELECTRICAL ENGINEER 


from the 
NAVAL POSTGRADUATE SCHOOL 
March 1971 





ABSTRACT 

An airborne direction-finding technique capable of 
handling multiple emitters was developed. An adaptive gate 
size was introduced in the track correlations based on co- 
variance relations of established tracks and observations. 
Track files were generated based upon the existence of 
various target parameters; i.e., frequency, pulse repetition 
frequency, pulse width and direction-of-arrival. To test 
the angular resolution capability of the filter, emitters 
in close proximity to each other with identical siete 
characteristics were used in the simulation. Target loca- 
tions are calculated in a carteSién coordinate system 
where the eee of the earth is taken into account, 
and with appropriate coordinate transformation computational 


Simplicity is preserved. 
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I. INTRODUCTION 


In airborne direction-finding (DF) a good procedure 
for locating the emitter given angle observations ‘only is 
of primary importance. This becomes very difficult ina 
multi-emitter environment. With sophisticated hardware and 
software, many of these emitters can be identified and sorted 
by their pulse meraminion frequency (PRF), and pulse width 
(PW). However, without accurate direction of arrival in- 
formation and methods for correlating these with successive 
past emissions, the problem is not solved. Central to all 
of this is the estimate of the location of the emitters. 

In the emitter locating portion of the study, it was 
assumed that all emitters have the same frequency, PRF, and 
pulse width, and thus, except for their angle of arrival, 
all correlate as the same target. This is really a very 
difficult problem. To get an idea just how difficult, 
consider a rather typical environment of six emitters in 
reasonable proximity, and an aircraft with a noisy antenna 
system making ten observations of each of these emitters. 
The possible number of locations with a straight-forward 
method of locating by calculating all possible intersections 
of bearing lines will give 1,620 locations computed by 


the formula 


mon (n=1) /2 


where m = number of emitters and n = number of observation 


points. The computer does not know whether there are Sj 
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emitters or 1,620. Without sophisticated association and 
correlation routines for correlating new observations with 
previous tracks and adaptive angle tracking and gating, 
the problem is virtually insolvable. 

Two methods of solving the location and correlation 
problem were investigated. The first uses a Kalman filter 
to compute predicted values of angle of arrival based on 
past data, and a correlation routine to correlate the 
present observations with these predicted values. If no 
correlation was possible, it initializes a new track. A 
correlated observation is then used to compute an optimal 
estimate of the present angle-of-arrival for that track. 
In the second method, all data fron a number of observation 
stations are processed simultaneously by first calculating 
all possible target locations based on intersections with 
the first bearing obtained. Then, the mean position of 
this cluster is found, declared a probable target location, 
and the observations having contributed to this location 
are correlated to this target. In both cases, the target 
location is found in an earth-centered coordinate system 
plane 


(where the x, axis is the polar axis and the x x 


ek 2 
is the equatorial plane of the earth, with the positive xy 
axls corresponding to 0° longitude and the positive X5 

axls corresponding to 90° longitude east). For each obser- 
vation, an observation plane is defined, passing through 


the earth's center and the bearing line at the observation 


point, such that in the noiseless case, this plane would 





also pass through the target point. The er plane is then 
defined by 


eel A SIGE OSG. a a 


aa ON, (1) 
where the coefficients as4 can be found from the observaticn 
station coordinates and the observation angle. The procedure 
for obtaining a target location is to take the equations of 
the observation planes, two at a time, with the constraint 
equation that the emitter be located at the earth's surface. 


Solving the equations of two planes below, gives the desired 


target location. 


a3 (2; 


j 't—- — 
xi +a = -a,,Cc (3) 


1 and X5 with 


X, = C. The constraint equation is accounted for by 


The equations were modified to solve for x 


SHarving for 


E = R/ (x"4 - aa ze 07) 1/2 


(4) 
where R = radius of the earth. The resulting solution 
Pope, 15 

a 


x, = E x! Sie ly 293. (5) 


By setting X, = C, the observation planes have been projected 
into lines on the plane X, = C, then in solving (2) and (3) 


the intersection of these lines in the X, = Cc plane is 
found. The parameter E is then utilized to project this 
point to the surface of the earth. With proper coordinate 


rotation, it is possible to have the projection plane (x,=c) 


earth surface tangent and aircraft centered. 
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IIT. COMPUTATIONAL PROCEDURE 


A. VARIOUS APPROACHES 

As 1s always the case, there are many ways to attack 
this problem. Possibly the best technique for evaluating 
the various alternatives is through a realistic statistical 
evaluation of a multi-emitter environment. But first a 
careful description of the two alternative techniques 
considered here, and the methods Be analyzing the tech- 
niques (i.e., the simulation scenario) is in coe Beles te 
a procedure for generating the target and receiver kinema- 
tics is necessary. Then, for the multi-emitter environment, 
a very soohisticated association and correlation routine 
for relating the observations with existing tracks is 
required. Lastly, the location of the emitters must be 
found, given either estimates of bearing angles or, alterna- 
tively, a large set of noisy bearing lines giving a set of 


noisy locations from which a location estimate is derived. 


B. AN ADAPTIVE BEARING ANGLE FILTER 

A Kalman filter[1] is used to obtain estimates and 
prediction of bearing angles for the next observation point 
to allow the correlation of observations with past data. 
Each observation and bearing line will have an unknown 
angular rate associated to it. This rate depends on the 


speed and heading of the aircraft, the bearing itself, 





and the sien range to the emitter for fixed or slowly 
moving targets. It is, however, approximately known; 

and this approximation is used in the initialization of 

the filter for each target. The initial uncertainity of 
angle and angular rate on track initialization is accounted 
for through the apriori covariance of error initialization 
in the filter. The filter also takes into account missing 
data points. That is, at any instant of look, at most 

only 50% of the targets might be emitting. Trackfile 
management, like initializing a new track, dropping-a track, 
or flagging a missing observation are done in the correla- 


£i0n section. 


The noisy observation z is assumed to be described 


by 


Za =ho +t vy, IO ey Zap eso, (6) 


ie is the true observation 


oe hy =a Oe ecanicd n= (6), Oy 


angle and angular rate, and Vy represents a white noise 
sequence. Since the observer is moving in time (the 

ch : : dif- 
subscript k denotes the k lniSiecuhe ened ame ) ec Cll. 


ference equation is used to approximate the system dynamics: 
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(en. = ek? ° ee Le 2, eels (7) 


No random forcing is included; this implies stationary 
targets and exact knowledge of the observer's position. 
The filter can be extended to take moving targets into 
account by augmenting the difference equation above, and 
using an appropriate oe ene Q matrix when initializing 
the recursive filter gain and covariance equations. The 
mean and variance of the measurement noise Sequence v, 1s 


given by: 


E(v,) = © for all k (8) 

and ree = RS a : (9) 
where Say = 0, 3#k 
J 1, 3=k 


The recursion equations for filtering are defined below, 


“A 


where © /5 denotes the estimate of the state Oy based upon 


measurement data (Zs, Zoreeecs Zen) 


5 
See kao 1 /K=1 (10) 
Sy = Be pn 7 Se KBR Red) | oS 


where g, are the filter gains. 


The recursive equations for the filter gain and 


covariance matrixes are defined bv: 
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T T =| 
Bye eet Be SB Pi fppea Be + R! (12) 
ee eye ew De 7-1 (13) 
P = 9 Pp oe + 14 
k+1/k ~ °k+1,k *k/k °k4+1,k * © (14) 


Qis a measure of random exitation in the process and 
for a stationary emitter is omitted. 

-No estimate of the apriori paiue of © is available. 
Therefore, the filter is initialized with the piers ee eisucee 


ment in the following fashion: 


8 eo (15) 


i 


CMe 
lI 


(w/R) sin(z,) : (16) 


In equation (16), v is the velocity of the observer, 
R the distance to the target, and Z4 the measured bearing 
angle. The only unknown in this equation is RK, but this can 
be approximated. A value of 150 nmi is assumed in the 
computer simulation. 


During simulation it was found that a simple ae k=] 
a 


a T 
matrix-of the form 6 = does not adequately 
kek | 0 1 


represent the system. This can be seen in Figs. 2 and 3 

where the mean errorsof bearing angle and angular rate 

are shown using the simple and an augmented te jeeil matrix 
a 


in a Monte Carlo simulation. The augmentation was found 
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Fic, 1, SIMULATION SCENARIO AT FIRST OBSERVATION POINT WITH 
THREE OF SIX TARGETS RANDOMLY EMITTING AT AN AVERAGE RANGE 
oF 150 NMI. 


INITIALIZATION EQUATION: 9 = (V/R) 5 1N2(9) 






Mean Error (degrees) 





Sa | a 


Fic. 2. MEAN ERROR OF BEARING ANGLE OBTAINED FROM A 
MoNTE CARLO SIMULATION ON THE KALMAN FILTER WITH A SINGLE 
TARGET AND PERFECT CORRELATION FOR A TOTAL OF 20 STATION POINTS, 


(J MEAN. ERROR OF THE MEASUREMENTS 


+ MEAN ERROR OF ANGLE ESTIMATES USING A NONLINEAR & MATRIX 
IN THE KALMAN FILTER PREDICTION EQUATION 


X MEAN ERROR OF ANGLE ESTIMATE USING A LINEAR © MATRIX IN 
THE KALMAN FILTER PREDICTION EQUATION 
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IG; 


4, MEAN ESTIMATE OF ANGULAR RATE AND MEAN RATE ERROR 


TRUE ANGULAR RATE 


MEAN ESTIMATE OF ANGULAR RATE AND MEAN RATE ERROR 
OBTAINED WHEN USING NONLINEAR 2& MATRIX 


MEAN ESTIMATE OF ANGULAR RATE AND MEAN RATE ERROR 
OBTAINED WHEN USING LINEAR 2 MATRIX 


iD 





by introducing an additional difference equation: 


“a a 


+ 


eee /k-1 * “k/k-1- (17) 


The addition of another variable, the second derivative of 
the bearing angle, does not greatly complicate the computa- 
tion; because for stationary targets (see Fig. 1), 8 can be 
expressed in terms of already defined variables. The 


augmented recursion equation for ace then becomes: 


aA aA aA 


(14+2*6 


Oy /k-1 ~ °k=1/k-1 *cot (6 


)T) . (ey) 


Fee c= 


This did improve the filtering results and eliminated the 
bias. 

The easiest and fastest way to compute an apriori 
target location using the results obtained by the Kalman 
filter is to compute the target location based on the first 
observation made and on the estimate of bearing angle 
obtained at a later observation site, after a number of 
observations have been made. However, the first observa- 
tion is an unfiltered measurement, and therefore may have 
quite a large error which will pearade the apriori target 
position seriously. To overcome this difficulty, the 
measurements are not only used in the Kalman filter to 
compute a current estimate, but also to obtain a smoothed 
estimate of the first observation. This smoothing is done 


by a technique proposed by H. E. Rauch[2]. The smoothed 
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a 


estimate of the initial observation after k measurements 


are taken is given by: 


aA An aA 


er ee Tk 8k Ac-2) - (19) 
D,» =D P ot pt (20 
L/k Wiel /k—teeok k-1 “Kk /k-1 20) 


Equation (20) shows that if the simple ¢@ matrix is used, 

Di /k is data independent and may b2 precomputed and stored. 
For this case, and using the initial values given in 
Appendix A, 


1 k=.) 


yee > 


0 il 


and since ony a smoothed value of the bearing angle and not 
the angular rate iS required, only a single scalar equation 
is needed to obtain a smoothed eStimate of the initial 
bearing angle. 

This leads to the following set of recursive filtering 


and smoothing equations in scalar form: 


a 
A a 


Oat = ORerseea + T8x-/KH1 

Oe /k-1 = 8-1 /kol (14276) 1 yoy cot (0, jay) 

OK /k o Ok Ag=1 P21) = Oe ea) (21) 
x = bes ke1 oye eee x pent) 

21 Ae - oye * Oye (1) (2,7 8 pny) 7 (K-1) Gy, (2) (24-8 pny) 
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e. PARAMETER CORRELATION 

The purpose of correlation is to assign observations 
to already established tracks if possible, or flag an obser- 
vation as a new target if this observation does not correlate 
to previously established tracks. To each observation (or 
track) five parameters are associated: 

Frequency 

PRF (Pulse Repetition Frequency) 

PW (Pulse Width) 

Correlation counter 

DOA (Direction of Arrival) 

For each of these paramters a gate is selected. The size 

of the gate depends on the expected accuracy of the measure- 
ment and the confidence associated with the track parameter 
in question. 

First, an observation is correlated to a track if the 
values of all five parameters are within their respective 
pategwmdiuhe This may correlate one observation with several 
tracks, or several observations with the same track. To 
eliminate this multiple correlation, one rejection rules 
and one control rule are used: 

1. A track which correlates with several observations 

rejects any observation held in common with another 
track if the common observation is the only obser- 
vation correlating with the other track. 

2. A track correlating with several observations, 


some of which are not held in common with other 


tracks, rejects observations held in common with 


other tracks. 
18 





3. When several observations correlate with one track, 
the closest observation is correlated to that track. 
Seeocomeseveral tracks Correlate with one observation, 
the observation is correlated with the closest 
track. 
5. If an observation is found uncorrelated, and at 
the same time there are cne or more tracks uncorre- 
lated to any observation, then this observation 
is correlated to the first uncorrelated track 
found within the gates. 
Biemgetes are chosen to be two Sioma of the target parame- 
Peed SUuLeMent MO1Se plus two sigma of the confidence inter- 
val associated with the track parameter estimate error 
which is a function of the number of observations already 
Sermroltated to that track. 
Rule 5 was included, because rules 3 and 4 can reject 
too many of the initial correlations. If, for instance, 
two observations are correlated to the same two tracks, 
it is possible that one of these observations--say, the 
Pest eaomcloser to both tracks than the second observation. 
Then rule 3 rejects the second observation from both tracks 
and leaves the first correlated to both tracks. Rule 4 
then selects the closest track as correlated to the 
first observation. This leaves unconnected the second 


Beservation and one track which was correlated initially. 
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D. TARGET LOCATION AND CORRELATION BY CLUSTER POINTS 

The intent of this location method is to find the most 
probable t:arget locations based on observations made at 
several stations, flag the observations leading to this 
location, and give a standard deviation to reflect the quality 
of the obtained observation. Fig. 4a shows four observaticn 
stations (Points I, II, III, IV) and possible bearings at 
these stations. A look at this Figure indicates quite 
Clearly that points X, Y, and Z are not only the most 
probable, but also the only logical choices for target 
locations. To find these locations by means of a digital 
computer, one observation from station I is selected 
(Fig. 4b), and crossing points with all bearings from all 
other stations are calculated. Then, the mean and standard 
deviation of these locations is found. Since in the close 
neighborhood of point X, three points are located; the mean 
position will be relatively close to this point, but the 
far away points will make the standard deviation large. 
Now all points outside a 1.lo interval from the mean are 
eliminated, and with the remaining points a new mean and 
Standard deviation is calculated. This procedure will elimi- 
nate erroneous bearing lines and will finally produce a 
mean very close to the point X with only the appropriate 
bearings from the other stations contributing to it. Now 
if the computed standard deviation is small enough to be 
Weccmtianea location confidence Limit, X 1s called a pro- 


bable target location. The observations Xr Xr X3r Xq are 
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flagged as correlated. These observations are deleted 

from further computations. Then the next observation from 
station I is selected and crossings with the remaining 
Bearing lines from all other stations are calculated 

(Fig. 4c). Here again, by computing a mean location and 
associated standard deviation, and rejecting distant cross- 
ings, the computed mean should converge to a location close 
to Y with a standard deviation smaller than the preassigned 
limit. Then for the last case, ideally only bearings lead- 


ing to point Z are left. 


Zen 





flagged as correlated. These observations are deleted 

from further computations. Then the next observation from 
station I is selected and crossings with the remaining 
Seeeang tines from all other stations are calculated 

(Fig. 4c). Here again, by computing a mean location and 
associated standard deviation, and rejecting distant cross- 
ings, the computed mean should converge to a location close 
to Y with a standard deviation smaller than the preassigned 
limit. Then for the last case, ideally only bearings lead- 


ing to point Z are left. 
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flagged as correlated. These observations are deleted 

from further computations. Then the next observation from 
station I is selected and crossings with the remaining 
bearing lines from all other stations are calculated 

(Fig. 4c). Here again, by computing a mean location and 
associated standard deviation, and rejecting distant cross- 
ings, the computed mean should converge to a location close 
to Y with a standard deviation smaller than the preassigned 
limit. Then for the last case, ideally only bearings lead- 


ing to point Z are left. 
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Peo kaooN LALO OF RESULTS 


A. TARGET LOCATION 

For the computer simulation, two target configurations 
were taken. In the first, three targets were placed 
approximately 18 nmi apart on a line parallel to the 
heading of the observer aircraft. In the second, six tar- 
gets were placed on a circle of approximately 50 nmi 
diameter, with the center of the circle approximately 
125 nmi from the path of the observer. As a final 
aggravation, the probability of detection of any target was 
reduced to 1/2; i.e., only three targets were emitting at 
any instant. A random sequence determined which three out 
of the six targets were observed at each station. 

In the simulation the "cluster" method was first tested, 
using target configuration one; Fig. 5 shows the results 
obtained. Only two out of the three targets were found. 
The position errors were less than ten miles. The third 
Eaueget could be recovered by taking a second look at the 
still uncorrelated observations. To get an estimate on 
the spread of the computed target locations that could 
be expected, a Monte Carlo simulation was done. The result 
is shown in Fig. 6; the greatest error was approximately 
50 nmi, and the center target was missed twice. The com- 
puted mean location from ten runs was on the far side for 


all three targets with an error of about 7 nmi. 
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Fic. 5, TARGET LOCATION BY “CLUSTER” METHOD, SINGLE RUN, 
THREE TARGETS ARE APPROXIMATELY 18 IM APART: 4 OBSERVATION 
POINTS ARE ALSO APPROXIMATELY I8NM APART. THE DISTANCE 
BETWEEN OBSERVATIONS AND TARGETS IS APPROXIMATELY 150 Ki, 


ALSO SHOWN ARE THE NOISY OBSERVATIONS MADE AT EACH STATION, 
+ TRUE TARGET LOCATION 


AW CALCULATED “CLUSTER” TARGET LOCATION 


X OBSERVATION POINT 
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Fic, 6, TEN RUNS OF TARGET LOCATION BY “CLUSTER” METHOD 


+ TRUE TARGET LOCATION 

A. COMPUTED TARGET LOCATION, INDIVIDUAL RUN 
Hi MEAN TARGET LOCATION FROMALL RUNS 

X OBSERVATION POINTS 
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Fic, 7, TARGET LOCATIONS OBTAINED FROM “CLUSTER” METHOD, 
SIX TARGETS ARE PLACED ON A CIRCLE OF 25 NM RADIUS. THE 
CENTER OF THE CIRCLE IS 125 NM FROM THE BASE LINE OF OBSER- 
VATION STATIONS, 


+ TRUE TARGET LOCATION 
4 COMPUTED TARGET LOCATION 
X LOCATION OF OBSERVATION POINTS 
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The results using the second target configuration show 
only one of the computed locations is reasonably close to 
the corressonding target, and that looking for cluster 
points is not very effective in a peenecer environment. 
The number of false locations becomes so high that the pos- 
Sibility of random clusters formed is too high to make the 
occurrence of clusters a good criterion for target location. 
Therefore, a new method mae developed using a Kalman filter, 
and a correlation procedure to obtain refined estimates of 
the observation angles from the measurements eae a 
The result is shown in Fig. 8; all positions except one are 
closer than 10 nmi to the true target position. To allow 
a comparison between the two methods, the same distance 
from the first to the last observation point was chosen. 

In Fig. 9 the tracking distance is doubled; this gave a 
considerable improvement. The largest error now was 4.3 nmi 
with 3 positions closer than 2 nmi to the true target 
positions. 

As a final test, the probability of detection was 
reduced to 1/2. Fig. 10 shows the results obtained; given 
on the graph are the initial Beene line (smoothed) where 
a new track was initialized; the estimates of bearing lines 
at the last station; the target location based on only 
these two bearings (0); the location based on all estimates 
available (1); the true target locations (+); and the 
observation points (X). The result reflect the missing 


observations. However, the computed target positions are 
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FIG. 8. TARGET LOCATIONS OBTAINED WHEN USING KALMAN FILTER 
ON OBSERVATION ANGLES. JHE BEARING LINES SHOWN ARE THE 
SMOOTHED ESTIMATES OBTAINED AT THE FIRST AND LAST OBSERVATION 
POINTS, 


t+ TRUE TARGET LOCATION 


1 MEAN TARGET LOCATION OBTAINED WHEN USING ALL BEARING 
ESTIMATES AVAILABLE, (2487+) Locations FOR EACH TARGET) 


Q TARGET LOCATION USING BEARING ESTIMATE OF FIRST AND LAST 
OBSERVATION POINT ONLY 


X OBSERVATION POINTS 
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Fic. 9. TARGET LOCATION USING KALMAN FILTER WITH 20 OBSER- 
VATION STATIONS. ALL TARGETS ARE OBSERVED AT ALL STATIONS. 


-+- TRUE TARGET LOCATION 


[} MEAN TARGET LOCATION USING ALL AVAILABLE BEARING 
ESI MAWES 


Q TARGET LOCATION USING FIRST AND LAST BEARING ESTIMATE 
ONLY 


X OBSERVATION POINTS 


Zo 





Degrees 
Latitude 
3.0 





Degrees 
Longitude 


Fic. 10. TARGET LOCATION USING KALMAN FILTER WITH 20 OBSER- 
VATION STATIONS. ONLY 3 TARGETS ARE OBSERVED AT ANY ONE 
STATION, THE GRAPH SHOWS THE ESTIMATE OF THE INITIAL OBSER- 
VATION AND THE LAST OBSERVATION OF EACH TRACK, 


+ TRUE TARGET LOCATION 
MEAN TARGET LOCATION BASED ON ALL AVAILABLE ESTIMATES 


OOD 


TARGET LOCATION BASED ON FIRST AND LAST BEARING ESTIMATE 
ONLY : 
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still close enough to the true positions to be useful as 


apriori target locations. 


B. TARGET CORRELATION 

A aomace of error arises in the target correlation. 
Two conflicting requirements make a compromise necessary. 
First, to reduce the possibility of correlating with wrong 
observations, it is desirable to select a small gate size, 
This, however, increases the possibility of missing observa-~ 
tions which should be correlated with a track but happen 
to be noisier than expected. This is undesirable, pecatee 
an uncorrelated observation will cause initialization of a 
new track very close to the already existing track. Now 
these two tracks will both be correlated to observations 
made from the same target, and these observations will be 
Spite by the two tracks--degrading the quality of both. 
So, false initialization of new tracks is even less desira- 
ble than allowing the possibility of a false correlation. 
tie meetiOunancewor the Correlation method was evaluated 
for several test runs. Tables I and II show two typical 
cases. The numbering of tracks and observations corresponds 
to the numbering of the targets shown in Fig. 7. 

The percentage of wrong correlations of 20% and 33.3% 
is quite high and degrades the filter performance seriously. 
Further studies should be done to find a more effective 
way of correlation. Even when noisefree observations 
were used, the noise introduced by the filter due to the 


approximation of the initial angular rate and the closeness 
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TABLE | 


CORRELATION RESULT FROM 20 OBSERVATION POINTS WITH ALL 
6 TARGETS OBSERVED AT EACH STATION. ALL TRACKS ARE 
INITIALIZED AT THE FIRST STATION. 

120 oBSERVATIONS 

24 WRONG CORRELATIONS ---™ 202 
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TABLE I] 


CORRELATION RESULT FROM 20 OBSERVATION POINTS WITH 3 
TARGETS OBSERVED AT EACH STATION. JHE NUMBER IN PAREN- 
THESIS GIVES THE STATION AT WHICH THE TRACK WAS INITIAL- 
IZED, - TRACK NUMBERS ARE ASSIGNED EQUAL TO THE TARGET 
NUMBER THAT HAS CONTRIBUTED THE MOST OBSERVATIONS TO 
THIS TRACK, 

60 OBSERVATIONS 

20 WRONG CORRELATIONS --*® 33,34 
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a 


of the targets considered was sufficient to cause correlation 
errors. This was particularly so in the low probability 

of detection case, when targets were emitting close together, 
and only one track was initialized. Then this track tends 

to pull in the observations from both targets, unless 


they are observed simultaneously at the same station. 


C. STATISTICAL EVALUATION 

Finally, a Monte Carlo simulation was performed using a 
Single target (andideal correlation) to verify that bias-free 
estimates are generated, and that the quality of aye esti- 
mates are properly represented by the analytical covariances 
obtained. The evaluation is done for the Kalman filter, 
Tearneana without the correction in the $(2,2) term. In 
both BAsaae the covariance sequences of the estimates of 
bearing angle and angular rate are reasonably close to 
the analytical values (Fig. 11, 12); but without the correc- 
tion term, a bias is built up in the estimate of bearing 
angle and the estimate of the angular rate. With the 
correction term, no bias is present, and the mean error 


of the estimate approaches zero (Fig. 2,3). 
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Angle Variance (degrees 





Fic. 11, VARIANCE OF BEARING ANGLE ESTIMATE OBTAINED FROM 
A Monte CARLO SIMULATION ON THE KALMAN FILTER WITH A SINGLE 
TARGET AND PERFECT CORRELATION. Fypst oBSERVATION IS AT 


Ee ZERO), 
(J) ANALYTICAL VARIANCE OF THE ANGLE ESTIMATE 


+ EXPERIMENTAL VARIANCE OF THE ANGLE ESTIMATE WHEN USING 
NONLINEAR & MATRIX IN KALMAN FILTER PREDICTION EQUATION 


X EXPERIMENTAL VARIANCE OF THE ANGLE ESTIMATE WHEN USING 
LINEAR ® MATRIX IN KALMAN FILTER PREDICTION EQUATION. 
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16. Time 

(min) 
Fig, 12. VARIANCE OF ESTIMATE OF ANGULAR RATE OBTAINED FROM 
A Monte CARLO SIMULATION ON THE KALMAN FILTER WITH A SINGLE 
TARGET AND PERFECT CORRELATION. 


[1 ANALYTICAL VARIANCE OF THE RATE ESTIMATE 


+ EXPERIMENTAL VARIANCE OF THE RATE ESTIMATE WHEN USING 
NONLINEAR ® MATRIX 


X EXPERIMENTAL VARIANCE OF THE RATE ESTIMATE WHEN USING 
LINEAR ® MATRIX | 
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IV. CONCLUSIONS 


The results from the computer simulation showed that 
a straight forward search for cluster points is not Pees 
tive in a high-density, multiple-emitter environment. This 
1s so becuase if m targets are emitting, and there are 
n observations, only 1/m of all possible locations are 
true locations--all others are wrong! This ratio of true 
to faulty locations can be reduced, and this was done in the 
computer simulation by first eliminating locations on the 
wrong side of the flight path, eliminating locations too 
distant ts. be considered meaningful, and by considering 
locations along selected bearings only. The high number of 
faulty locations still remaining is sufficient to generate 
erroneous clusters to which the process may converge and 
lead to false locations. 

When a Kalman filter was used to process the measure- 
ments, the critical point was the correlation procedure 
which selects the measurements and assigns them to a track. 
If this selection was done perfectly, the filter would have 
only the proper data as its input and would produce an opti- 
mal estimate of the systems state. As correlation errors 
occur, wrong data is fed into the filter and performance 
deteriorates. This indicates that the overall system is 
only as good as the correlation. For this reason, the 
correlation procedure is quite elaborate; adaptive gating 


was introduced and five rules are used to handle the 
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multiple correlation problem as efficiently as possible. 
However, it must be pointed out that the two requirements 

for a good correlation method--correlation and discrimination-- 
in a sense contradict each other and cannot both be maximized 
at the same time. In the initial correlation, a greater 
weight 1S put on the correlation; and the 4 rejection rules 
provide discrimination, but errors do occur during this 
rejection process, and the high error percentage obtained 
Suggests that a better correlation method could be found 

which then, in turn, would improve the filter performance 


and the quality of locations obtained. 
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APPENDIX A 


‘KALMAN FILTER BLOCK DIAGRAM AND INITIALIZATION 


Recursive gain and covariance equations for filtering: 


A a -l 
Sx = Py yx hy xP Ksx-Ext Rx! 
PK/K e/a Dey /K-1 
P = 9 Ee oe + Q 
K+1/K  ~“K+1,K K/K K+1,K K 
Initial values: 
10° 0 1 Tt 
WAG ee 
0 ip t 
0 1 
h, = [1 0] 
Ri = l 
0 0 
Ox ~ for stationary targets. 
0 0 
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Recursive 


D 


STK 


Initial values: 


D 


e 


lek 


iy 


OF 


gain and covariance equations for smoothing: 


T ae 


ea em Ko K/K=1 


Z _ T 
ee ee eo Re 1 7K 


Initial state. values: 


=| 
iy 
“sin z 
R II, 
Ve=eve locity Of alrcrafrt 
R = 150 NM assumed initial distance to target 
Note: In the Kalman filter equation a non-linear 6 matrix 


is used;when calculating gains and covariances, the linear 


approximation of the $ matrix is used. 
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Kalman filter gain and covariance equations in scalar form. 


dy (1) 


Jy, (2) = 


Peni 


PK+L/K 


PRAL/K 


i PK/K-1 tI IPR py 


Rete ee 


(1,1) = g,(1)R + 9,(2) [2R-P 


ey ae 
(1,1) + R} 


ey 102) + Pigg, (252) 


(172) = Gy (2) [R-Py ppg (102) 1 + Py py (252) 


(2,2) = Py py 3 (212) - Py yyy (11209 (2) 
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Fig, 13. BLockK DIAGRAM OF KALMAN FILTER 
AND SMOOTHING SOLUTION 
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COMPUTER RESULTS FOR SIX TARGETS AND 50% DETECTION 
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multiple emitters was developed. An adaptive gate size was introduced 
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